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ABSTRACT 

We report the detection of CO (3— >2) emission from the submillimeter-selected luminous galaxy 
SMM J14011+0252. The optical counterpart of the submillimcter source has been identified as a merger 
system with spectral characteristics consistent with a starburst at z = 2.565. The CO emission confirms 
the optical identification of the submillimeter source and implies a molecular gas mass of 5xlO lo /if 5 2 M , 
after correcting for a lensing amplification factor of 2.75. The large molecular gas mass and the radio 
emission are consistent with the starburst interpretation of the source. These results are similar to 
those found for SMMJ02399-0136, which was the first submillimeter selected CO source found at high 
redshift. The CO detections of these two high-redshift submillimeter galaxies suggest the presence of 
massive reservoirs of molecular gas which is consistent with the inferred high rates of star-formation 
(10 3 Mq yr _1 ). These two systems appear to be associated with merger events which may evolve into 
present day luminous elliptical galaxies. 

Subject headings: early universe — galaxies: active — galaxies: evolution — galaxies: individual 
(SMM J1401 1+0252) — galaxies: starburst 



1. INTRODUCTION 



2. OBSERVATIONS 



Deep surveys of the submillimeter sky using the Submil- 
limeter Common-User Bolometer Array (SCUBA) on the 
James Clerk Maxwell Telescope have uncovered a popula- 
tion of distant dust-rich galaxies (Smail, Ivison & Blain 
1997; Blain et al. 1998; Barger et al. 1998; Hughes et al. 
1998; Eales et al. 1998). Based on their optical colors, the 
majority of these systems are thought to lie at redshifts of 
z ~ 1-5 (Smail et al. 1998b). However, only a few of these 
systems currently have accurate redshifts. The first well- 
studied system SMM J02399-0136, hereafter SMM J02399, 
at z = 2.8 was shown to contain both an AGN (Ivison et 
al. 1998) and a massive reservoir of molecular gas thought 
to be fueling a starburst (Frayer et al. 1998). 

In order to test the generality of the properties of 
SMM J02399, we have initiated a search for CO emission 
in additional sub-mm sources. In this Letter we report the 
results for the submillimeter source SMM J14011+0252, 
hereafter SMM J14011. This sub-mm galaxy was discov- 
ered during a survey through rich, lensing clusters (Smail 
et al. 1998b). Follow-up optical spectroscopy places the 
sub-mm source at a redshift of approximately z — 2.55 
(Barger et al. 1999; Ivison et al. 1999). SMMJ14011 ap- 
pears to be a merger/interacting system which has optical 
and radio properties consistent with an ultraluminous star- 
burst (Ivison et al. 1999), albeit the presence of a heavily 
obscured AGN cannot be ruled out. The amplification fac- 
tor for SMM J14011 is 2.75 ± 0.25, where the uncertainty 
comes from the allowed range of detailed mass models of 
the foreground cluster. 



SMM J14011 was observed using the Owens Valley Mil- 
limeter Array between 1998 October and December. A to- 
tal of 41 hours of usable integration time on source was ob- 
tained in two configurations of six 10.4m telescopes. The 
phase center for the CO observations was the position of 
the brightest optical component of the SMM J 140 11 sys- 
tem (Jl); a(J2000)=14 h 01 m 04?97; 5(J2000)=+02°52'24'/6 
(Ivison et al. 1999). The CO (3^2) line was observed us- 
ing a digital correlator configured with 112 x 4 MHz chan- 
nels centered on 96.9975 GHz in the lower side-band, cor- 
responding to CO (3— >2) emission at the Ha redshift of 
z — 2.565 (Ivison et al. 1999). Typical single sideband 
system temperatures were approximately 300-350 K, cor- 
rected for telescope losses and the atmosphere. In addition 
to the CO line data, we recorded the 3mm continuum data 
with a 1 GHz bandwidth for both the upper (line-free, cen- 
tered on 99.9975 GHz) and lower sidebands. The nearby 
quasar 1413+135 was observed every 25 minutes for gain 
and phase calibration. Absolute flux calibration was deter- 
mined from observations of Uranus, Neptune, and 3C 273. 
The absolute calibration uncertainty for the data is ap- 
proximately 15%. 



3. RESULTS 



Figure 1 shows the CO(3-^2) spectrum for SMM J14011. 
The CO line is detected at the redshift of the narrow Ha 
line of SMMJ14011 (Ivison et al. 1999). We achieved a 
7cr detection for the peak in the integrated CO (3-^2) map 
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(Fig. 2). The CO position (Table 1) is consistent with the 
optical position of SMMJ14011 within the uncertainties 
of the data sets. The CO map shows marginal evidence 
for emission extended toward the south, but observations 
at higher resolution are required to test this possibility. 
The upper limit of the continuum emission is S„(3mm) 
< 0.6 mJy (ler). This is insufficient to detect the thermal 
dust emission discovered by SCUBA, assuming S v oc f 3 5 . 

Both the CO and Ha lines of SMM J14011 are redshifted 
from the UV lines by about 1100 ± 700 km s" 1 (Ivison et 
al. 1999). Systematic blue-ward offsets of UV lines are 
not unusual for starbursts and have been attributed to 
outflows with dust obscuration at systemic and redshifted 
velocities (Mirabel & Sanders 1988; Gonzalez Delgado et 
al. 1998; Heckman et al. 1998). In an extreme example, 
the z = 3.9 quasar APM 08279+5255 shows a 2500 km s" 1 
blue-ward offset of the UV lines from the systemic CO red- 
shift (Downcs et al. 1998). 

The observed CO(3^2) line flux is S(CO) = 2.4 ± 
0.3 Jykm s _1 . No adjustment has been made to account 
for the continuum level since it is negligible. The observed 
CO (3— >2) line flux implies an intrinsic CO line luminosity 8 
of L'(CO) = 1.2 x lO^hfj; 2 Kkm s" 1 pc 2 (sec formulae in 
Solomon, Downes, & Radford 1992). The CO luminosity 
is related to the mass of molecular gas (including He) by 
M(H2)/L'(CO) = a. The value for a is expected to be be- 
tween a ~ lM©(Kkm s _1 pc 2 ) -1 (Solomon et al. 1997) 
and the Galactic value of a ~ 5M (Kkm s" 1 pc 2 )" 1 
(e.g., Sanders, Scoville, & Soifer 1991). We adopt a 
value of a = 4 M Q ( K km s _1 pc 2 )" 1 , which is consis- 
tent with that estimated for Arp 220 (Scoville, Yun, & 
Bryant 1997a) after correcting for the line brightness ratio 
of T 6 (CO[3 -> 2])/T b (CO[l -> 0]) ~ 0.6 typically observed 
in starbursts (Devereux et al. 1994). The inferred molec- 
ular gas mass of SMM J14011 is 5 x lO lo /if 5 2 M , which is 
consistent with that of the most massive low-redshift ul- 
traluminous infrared galaxies (ULIGs; Sanders & Mirabel 
1996). 

The implied gas-to-dust ratio for SMMJ14011 is 
M(H 2 )/Mdust = 200-900, assuming a reasonable range of 
possible dust temperatures (30-70 K). These gas-to-dust 
ratios are similar to those seen in spiral galaxies (Devereux 
& Young 1990), local ULIGs (Sanders et al. 1991), and 
other high-redshift CO sources. The similarities of the de- 
rived gas-to-dust ratios at high-redshift reflect the small 
scatter (only a factor of 2-3) in their L'(CO)/L^(submm) 
luminosity ratios. These results are somewhat remarkable 
considering that some of the sources are more likely domi- 
nated by AGNs than others and considering the possibility 
of different dust temperatures. 

The observed line width for SMMJ14011 
(FWHM=200km s" 1 ) is much narrower than that ob- 
served for SMMJ02399 (700km s" 1 ), but is within the 
range of values typically observed for ULIGs (200- 
400km s -1 ; Solomon et al. 1997). The narrow line - width 
may indicate that we are viewing the system somewhat 
face-on. By using the observed line width and the upper 
limit on the CO source size of 9 < 7", we can constrain the 
total dynamical mass contained within the CO emission 
regions. The angular size limit corresponds to a maxi- 

8 We adopt q = 0.5, H = 75 /175 km s _1 Mpc -1 , and a lensing a: 



mum linear diameter of D < 36/i 75 kpc at z — 2.565, or 
D < 13h^ kpc after correcting for lensing. The dynami- 
cal mass is M dyn ~ R(AV/[2 sin(i)]) 2 /G, where AU is the 
observed FWHM line width and i is the inclination. We 
find Mdyn < 1-5 x 10 10 sin _2 (z) h^ 5 Mq which is consistent 
with the derived gas mass, provided i <J 35°. 

4. DISCUSSION 

The intrinsic far-infrared luminosity of SMM J 140 11 is 
approximately 3 x 10 12 h^ L Q . If its luminosity is powered 
by star-formation, as suggested by the optical data (Ivi- 
son et al. 1999), the implied star-formation rate of mas- 
sive stars (M > 5M Q ) is SOOMoyr" 1 (Condon 1992). 
Including the presence of low mass stars, the total star 
formation rate would be of order 10 3 MQyr _1 . The CO 
data is consistent with this view by showing the pres- 
ence of enough gaseous material to support a high level 
of star formation. The large molecular gas mass and the 
high star-formation rate suggest a starburst time scale of 
order 5 x 10 yr. In addition, the far-infrared-to-radio 
flux ratio for SMMJ14011 is the same as that found for 
nearby starburst galaxies (Condon, Frayer, & Broderick 
1991). Alternatively, if SMMJ14011 is powered by an 
AGN, the AGN could not be a powerful radio source 
(< 3 x lO 30 ^ 2 ergs" 1 Hz" 1 ). 

The success of the first two searches for CO from the 
high-redshift submillimeter-selected galaxies is encourag- 
ing. Unlike SMM J02399, whose UV/optical spectral prop- 
erties show significant AGN activity, there is currently no 
evidence for an AGN in SMM J14011. However, the pres- 
ence of a dust-enshrouded AGN cannot be completely 
ruled out. Even though the UV/optical characteristics 
of SMMJ14011 and SMMJ02399 appear different, their 
radio, submillimeter, and CO properties are fairly simi- 
lar and suggest that star formation is important in the 
far-infrared emission of both sources. SMM J 140 11 and 
SMM J02399 share many of the same properties of the lo- 
cal population of ULIGs, such as (1) high infrared lumi- 
nosities, (2) associated with mergers, (3) large reservoirs 
of molecular gas, and (4) comparable CO line widths. The 
most important distinction is that the high-redshift sub- 
mm galaxies are about a factor of 10 2 -10 3 times more 
numerous per unit co-moving volume than the low red- 
shift ULIGs (Kim & Sanders 1998). The sub-mm galaxies 
may also have slightly larger amounts of molecular gas (by 
factors of 2-3) than typical ULIGs. 

Figure 3 shows the CO luminosities plotted as a func- 
tion of redshift. When lensing corrections are taken into 
account, the high-redshift sources have similar CO lumi- 
nosities to the most luminous low redshift ULIGs. These 
consistencies could suggest roughly similar gas masses in 
the progenitor systems or may indicate the importance of 
self-regulating mechanisms such as galactic winds (Heck- 
man, Armus, & Milcy 1990). For comparison, theoret- 
ical evolutionary curves are also presented in Figure 3. 
We show three models based on numerical calculations 
of chemical evolution (Frayer & Brown 1997). The mod- 
els are designed to represent a spiral galaxy, an ellipti- 
cal galaxy, and a merger event. All three models start 
as purely gaseous systems with zero metallicity. The spi- 
ral galaxy model uses an infall model (Run 11 of Frayer 

cation factor of 2.75 throughout this Letter. 
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& Brown 1997) which starts at z = 5 and reaches a fi- 
nal mass similar to that of the current Milky Way disk 
(6 x 10 10 Mq). The CO luminosity of a disk with the mass 
of the Milky Way is never expected to reach the high CO 
luminosities (;> 10 10 K km s~ 1 pc 2 ) observed in the high- 
redshift systems. For the elliptical galaxy model, we as- 
sume a closed-box model (Run 8 of Frayer & Brown 1997) 
with a mass of 2 x 10 11 Mq which begins its evolution at 
z = 5. Although galactic systems are not expected to 
evolve in such a simple manner, the closed-box model is 
able to match the high CO luminosities detected at high 
redshift and the low CO luminosities (<; 10 7 K km s _1 pc 2 ) 
observed for local elliptical galaxies (Wiklind, Combes, & 
Hcnkel 1995). The third model in Figure 3 is a more real- 
istic model for SMM J14011 and SMM J02399. It assumes 
a starburst resulting from the merger of two 2 x 10 11 Mq 
gaseous systems at z = 3. This starburst model uses the 
parameters of Run 8 in Frayer & Brown (1997), except 
for a factor of 10 increase in the star formation efficiency 
and assumes a formation time scale of r/ = 3 x 10 8 yr 
for the merger (Mihos & Hernquist 1996). These theo- 
retical calculations only show the approximate evolution 
expected for different types of galactic systems and depend 
strongly on the parameters governing star formation, the 
initial mass function, infall, and cosmology. 

An important unanswered question is what will 
SMMJ14011 evolve into at the present epoch. By using 
the observed near-infrared K-band flux density (Ivison et 
al. 1999) and assuming the burst model from Bruzual & 
Chariot (1993), SMM J14011 has already formed of order 
10 10 Mq of stars. There is enough molecular gas present 
to form an additional L* galaxy (adopting 5 x 10 10 Mq 
per L*), suggesting that SMMJ14011 will evolve into a 
~ 1-2 /iy 5 2 L* galaxy. Given that the low-redshift ULIGs 
may eventually evolve into elliptical galaxies (Kormendy & 
Sanders 1992; Mihos & Hernquist 1996), the high-redshift 
sub-mm population could represent luminous elliptical 
galaxies in their formative phases (Smail et al. 1997; Eales 
et al. 1998). The consistency of this hypothesis can be 
checked by comparing the number density of the sub-mm 
population represented by SMM J14011 and SMM J02399 



with that of luminous ellipticals at low redshift. Within an 
effective source plane area of 25 square arcminutes, Smail 
et al. (1998b) have detected 9 distant, non-cluster, sub- 
mm sources above the 4cr noise level. SMM J 140 11 and 
SMMJ02399 represent about 20% of this sample. The 
corresponding number density of such sources per unit co- 
moving volume within the redshift range of z — 1-5 is 
approximately 10~ 4 h^ 5 Mpc~ 3 . The number density of lu- 
minous elliptical galaxies found at low redshift is roughly 
similar (Loveday et al. 1992). Although no definitive con- 
clusions can be made due to the small sample size, the 
current data are at least consistent with the sub-mm pop- 
ulation evolving into the luminous elliptical galaxies at the 
present epoch. 

5. CONCLUSIONS 

We report the detection of CO in a second ultraluminous 
galaxy, SMM J14011, selected from a submillimeter survey 
of the distant Universe (Smail et al. 1998b). The CO emis- 
sion is coincident in position and redshift with the optical 
counterpart of SMM J 140 11 and indicates the presence of 
approximately 5 x 10 10 /iy 5 2 Mq of molecular gas. Both 
SMMJ14011 and the first high-redshift sub-mm galaxy, 
SMM J02399, show CO emission and are thought to be as- 
sociated with a major burst of star formation. The massive 
molecular gas reservoirs of SMM J14011 and SMM J02399 
and their co-moving number densities are consistent with 
these sources evolving into present day luminous ellipti- 
cal galaxies. Further observations are needed to constrain 
the redshift distribution and the molecular gas masses of 
sub-mm population in order to test the generality of these 
early results. 
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Fig. 1. — The CO(3^2) spectrum for SMM J14011 observed with the OVRO Millimeter Array. The data have been smoothed to 24 MHz 
(74 km s ) and the channel separation is 12 MHz. The la rms error of each channel is shown in the lower left. 

Fig. 2. — The integrated CO(3— >2) map averaged over 104 MHz (322 km s" 1 ) overlaid on an optical /-band image taken at the Palomar 
5 m (Smail ct al. 1998a). The positional offsets arc relative to the optical position, and the cross marks the position of the SCUBA source 
(Ivison et al. 1999). The la rms error is 0.32 Jy km s — 1 /beam, and the contour levels arc lax (—3,— 2,2,3,4,5,6,7). The synthesized beam size 
for the observations is shown in the lower left (6'.'5 X 4'.'3, PA = —26°). 

Fig. 3. — The L'(CO) luminosity as a function of redshift. The molecular mass scale (at right) assumes a = 4M Q (Kkm s-ipc 2 )" 1 . 
Besides SMM J14011 and SMM J02399 (boxes), the other high-redshift CO sources (large dots) are IRAS F 10214+4724 at z = 2.3 (Brown & 
Vanden Bout 1991; Solomon et al. 1992), 53W002 at z = 2.4 (Scoville et al. 1997b), the z = 2.6 Cloverleaf quasar H 1413+117 (Barvainis ct al. 
1994), MG 0414+0534 at z = 2.6 which has an unknown lensing magnification factor signified by the downward arrow (Barvainis et al. 1998), 
APM 08279+5255 at z = 3.9 (Downes ct al. 1998), BRI 1335-0417 at z = 4.4 (Guilloteau ct al. 1997), and BR 1202-0725 at z = 4.7 (Ohta et 
al. 1996; Omont ct al. 1996). For known lensed sources, the observed and intrinsic values are connected by a dashed line. In addition, we plot 
(small dots) the low redshift luminous IRAS sources from Sanders et al. (1991) and additional ultraluminous IRAS sources from Solomon et 
al. (1997). The Milky Way Galaxy is shown as a diamond (Solomon & Rivolo 1989). Theoretical evolutionary curves are shown as dotted 
lines (see text): (1) galactic disk model, (2) closed-box elliptical galaxy model, and (3) starburst merger model. 



Table 1 
CO Observational Results 



Parameter 


Value 


a(J2000) 


14 h 01 m 04?92 ± 0?03 


<5(J2000) 


+02°52'25"6 ± Of! 6 


z(CO) 


2.5653 ±0.0003 


Linewidth (FWHM) 


200 ±40 km s _1 


5(CO) a 


2.4±0.3Jykm s _1 


L(CO) b 


1.6 x lO 7 /if 5 2 L 


L'(CO) b 


1.2 x 10 10 /if 5 2 Kkm s^pc 2 


M(H 2 ) b < c 


5 x lO lo /i 75 2 M 



a Observed CO (3^2) line flux. 

b Intrinsic value assuming a lensing amplification factor of 2.75, q = 
0.5, and H Q = 75 h 75 km s _1 Mpc" 1 . 

c Estimatcd using a = 4M (Kkm s^pc 2 )" 1 . 
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